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ABSTRACT 

A fully coupled mathematical model describing the interactions between a vibrating thin 
cylindrical shell and an enclosed acoustic field is presented. Because the model will ultimately 
be used in control applications involving piezoceramic. actuators, the loads and material contri- 
butions resulting from piezoceramic patches bonded to the shell are included in the discussion. 
Theoretical and computational issues lead to the consideration of a weak form of the model- 
ing set of partial differential equations (PDE’s) and through the use a semigroup formulation, 
well-posedness results for the system model are obtained. 
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1 Introduction 


The active control of noise generated by structural vibrations has been studied for several 
years in various applications. One recent motivation leading to an intense study of prob- 
lems involving the reduction of structure-borne noise has resulted from the development 
of a new class of turboprop and turbofan engines. These engines, although very fuel ef- 
ficient, are also very noisy. Specifically, the low frequency, high magnitude exterior noise 
fields produced by these new engines cause vibrations in the fuselage which, due to struc- 
tural acoustic coupling between the fuselage and the interior acoustic field, lead to un- 
acceptably high cabin noise levels. The problem is exacerbated by the increased use of 
lightweight, composite materials in the cabin walls. In order to make the use of the en- 
gines feasible in commercial aircraft, a great deal of research has been aimed at developing 
control techniques to reduce this unwanted interior noise. 

Control techniques for structural acoustics problems of this type have been studied from a 
variety of perspectives [2, 7, 12, 13, 14, 15, 16, 20, 21, 25], one of which is through the use of 
piezoceramic patches which are bonded to the enclosing structure [2, 7, 14, 16]. As detailed 
in [8, 9], the patches create bending moments and in-plane strains when a voltage is applied, 
and through these actions, can be used to alter the structural dynamics in a manner which 
ultimately reduces interior noise. 

In this paper, we develop a fully coupled mathematical model, as well as well-posedness 
results for this model, which can be used when applying parameter estimation and PDE-based 
control strategies to problems involving the use of piezoceramic actuators. This model extends 
the 2-D results in [2, 7] and the 3-D results in [10] to the cylindrical domain of interest (the 
3-D domain in [10] consisted of a hard-walled cavity with a thin plate at one end). While 
ultimately motivated by the cabin noise problem mentioned above, the model is designed to 
be consistent with an experimental apparatus being used at the Acoustic Division, NASA 
Langley Research Center. This setup consists of a thin-walled hollow cylindrical shell which 
is supported by rigid caps at the ends. Piezoceramic patches bonded to the walls of the shell 
are then used to control interior noise which has been generated by the vibrations of the walls. 

This mathematical model differs from previously used 3-D shell/acoustic models in that 
it fully incorporates the backpressure and momentum conditions which couple the structural 
dynamics with the interior acoustic response, thus yielding a time-dependent system of partial 
differential equations which describe the system dynamics. Because the acoustic/structure 
interactions are fully incorporated, this model is useful for control strategies which utilize the 
natural “feedback” loop due to the coupling between the acoustic fields and the structural 
vibrations. Moreover, by incorporating the piezoceramic patch/shell interaction results of [8], 
the model can be employed in devising control techniques utilizing piezoceramic actuators. 

After the development of the mathematical model, well-posedness results are presented. 
These are important not only for determining the existence and continuous dependence of 
solutions, but also in providing an initial framework that can be used when determining 
suitable approximation schemes. By carefully noting the Hilbert spaces containing the state 
variables and test functions, appropriate choices can be made for the approximating subspaces 
which contain the basis functions used in finite element or spectral approximations. This is 
important not only when considering forward simulations but also when considering theoretical 
issues concerning the parameter estimation and control problems. 
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Section 2 of this presentation contains a description of the components of this system 
along with the strong form of the equations of motion for the system. The disadvantages of 
the strong form are discussed in Section 3 and a weak or variational form of the system equa- 
tions is presented. An abstract framework amenable to forward simulations, the estimation 
of physical parameters, and the implementation of PDE-based feedback control strategies is 
also developed in this section. In Section 4, the abstract model described in the third section 
is shown to be well-posed. This is accomplished using “extrapolation space” ideas and argu- 
ments similar to those presented in [5, 6, 17]. Having obtained the existence and continuous 
dependence of solutions for the model, approximation and LQR optimal control techniques 
similar to those discussed in [2, I, 10] can be applied to the problem in an effort to suitably 
reduce interior pressure levels. 


2 Strong Form of the System Equations 

As motivated by the experimental setup described in the Introduction, the structural acoustics 
problem under consideration is assumed to consist of a cylindrical acoustic cavity f l(t) which 
is enclosed by a thin cylindrical shell (see Figure 1). Hard wall conditions are taken at the 
ends F of the cavity in order to model the rigid end caps used in the experimental apparatus. 
These end caps also dictate the use of clamped boundary conditions when approximating the 
shell dynamics. 

As outlined in the discussion in the previous section, the dynamics of the coupled system 
are composed of an acoustic response, shell dynamics, piezoceramic patch/shell interactions, 
and the coupling acoustic/structure interactions. Each of these components is briefly described 
and the results are summarized in a coupled system of PDE’s which describe the system 
dynamics. 



Figure 1. Cylindrical acoustic cavity f 1(f) with hard end caps T. 

2.1. Acoustic response 

The acoustic wave motion inside the cavity can be described either in terms of a velocity 
potential 4> or the acoustic pressure p (the two are related through the relationship p = 
Pj(f> t where p f is the equilibrium density of the atmosphere); motivated by control theoretic 


2 


considerations and to simplify the presentation which follows, we choose the former as the 
second-order state variable. The acoustic dynamics inside the cavity are then modeled by the 
undamped wave equation 

<f) tt = c 2 A<f> , (r,$,x) € Q(t) ,t > 0, 

V<^ • h = 0 , ( 7 -, 0, x) 6 T , t, > 0 

where c is the speed of sound in the cavity, n is the outward axial unit normal to the cavity 
end caps, and the Laplacian in cylindrical coordinates is given by 

A«s = ^ + I^ + i^: + ^ 

dr 2 r dr r 2 dd 2 dx 2 

Again, T denotes the circular caps at the ends of the cavity Q(t) (the cavity is variable in 
time due to the fact that the enclosing shell is vibrating), and x is the axial direction in the 
cavity. Finally, we note that air damping inside the cavity was omitted due to the relatively 
small dimensions of the experimental cavity, and could readily be incorporated for significantly 
larger cavities. 

2.2. Shell dynamics 

In order to specify the equations of motion of the bounding structure, we consider a thin 
cylindrical shell having length l, thickness h and radius R with axial direction taken along 
the x-axis (see Figure 2). The displacements of the middle surface in the axial, tangential and 
radial directions are taken to be and w , respectively. Furthermore, the shell is assumed to 
have density p, Young’s modulus £?, Poisson ration //, and damping coefficient Co- Finally, wc 
assume that the shell’s length is relatively short in relation to it radius and hence the Donnell- 
Mushtari equations can be used when approximating its motion [18]. This last assumption 
is made purely for ease of presentation and higher-order theories such as the Byrne-Fliigge- 
Lur’ye model can be substituted for that of Donnell and Mushtari as warranted by the shell 
dimensions. 
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As presented in [8, 19], the Donnell-Mushtari equations for a thin cylindrical shell are given 


d 2 u dN x dN ex „d{N x ) pe 
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Here M x , Mg, N x and Ng are internal moments and force resultants, ( M x ) pe , (Mg) pe , (N x ) pe 
and (N$) pe are the resultants for the loads generated by the piezoceramic patches when a 
voltage is applied, and (j n is the external, normal load on the shell. For a patch with uniform 
thickness and bounding values X\,x 2 ,B\ and 0 2 as depicted in Figure 2, the presence of the 
indicator function 

1 , x < (*i + x 2 )/2 

S\, 2 {x) = < 0 , x = (x l +x 2 )/2 (2) 

-1 , x > (X! + x 2 )/2 

derives from the fact that the forces generated by the patch in the ar-direction are antisym- 
metric (equal in magnitude but opposite in sign) about the line x = (aq 4- a; 2 )/2. The same 
holds true for the forces in the 0-direction and Si, 2 (0) is defined in an analogous manner. 

From the discussion in [9], internal damping can be incorporated in the shell equations by 
assuming a constitutive law which posits that stress is proportional to a linear combination of 
strain and strain rate. This yields a a Kelvin- Voigt type of damping in the shell with internal 
resultants for shell regions not covered by patches given by 
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Similar expressions (e.g., see [9]) can be derived for those regions of the structure in which 
piezoceramic. patches are bonded to the shell. In those areas, the internal force and mo- 
ment resultants contain contributions not only from the shell but also from the internal 
stresses in the piezoceramic patches. The contributions resulting from the external loads 
(M x ) pe ,(Me) pe i(N x ) pe , (No) pe and q n will be discussed further in the next two subsections. 

Finally, we note that the boundary conditions for the shell must be consistent with the 
rigid end cap conditions on the experimental apparatus. Depending on the exact nature of 
the end caps, one may enforce either the clamped edge condition 


d w 

u = v = w = — — : 

ox 


0 


or the simply supported edge condition 


u = v = w = M x = 0 

at the ends of the shell. For definiteness, we will employ the first while simply noting that 
the latter often adequately models the conditions being maintained at the edges of the shell. 


2.3. Piezoceramic patch/cylindrical shell interactions 

In order to specify the moment and force resultants which result from the activation of the 
piezoceramic patches, we consider patches which are bonded to the inner and outer surfaces 
of the shell in a manner such that their edges are parallel to lines of constant x and 9 as shown 
in Figure 2. When it is necessary to differentiate between the two patches, the outer will be 
denoted with a subscript pe i with a subscript pe 2 being used to denote the inner patch in 
each pair. We point out that in order to keep this presentation compact, we have assumed 
that both patches are active, and we refer the reader to [8] for the case when only one patch 
is bonded to the shell. 

As discussed in [8], the total line moments and forces for a single patch pair with bounding 
values Xi,Z 2 , 0 i and d 2 are 

(M x ) pe = [(M$) pej + (M,) pe J [H,(x) - H 2 {x)\ [//,(£) - tf 2 (0)] 

(Mg) pe = \{M x ) pei + ( M x ) P e 2 ] [//i(z) — H 2 (x)\ [H\(6) — H 2 (9)\ 

( 4 ) 

(N x )pe — [(A x ) pe , + (Afr)y,e 2 ] [Hl(x) — H- 2 [x)] [//](<?) — H 2 (9)] Si >2 (x)S\ i2 (d) 

(A r e)pe — [(A^)pe, + [/fl (x) — H 2 (x )] [H\(9) — H 2 {6)] S\ i2 (x)Si i2 (6) . 

Here H is the Heaviside function with H t (x) = H(x - aq) ,i = 1,2, and S it2 is the indicator 
function defined in (2) with similar definitions being used in 9. The individual components 



are given by 


(Mx)pt i 



( 1 ”1“ ^pe) f-pej 


( Mx)pe2 


E. 


pe 


l - // 2 

pe 





( 1 ~b ^pe) ^pe2 


(A^)pei 



( 1 H - ^pe) ^pej 


(Mo)pe 2 
(Nx)pe i 


E t 


pe 


pe 


5N(j + r) 2 -^ 


( 1 H“ ^pe) ^ 


P?2 


(5) 




pp 


1 - // 


pe 


r + 7i5 4 (l + r ) -' ,2 


( 1 “b ^pe) 


PC I 


(^x)pc 2 



( 1 “b ^pe) e pP2 


(A^o)pe, — — . _ P 2 P ( 1 + »V) e Pfi 

* ^PP 

('' , «k = -TTTJ’(‘+v)' w 

* ^pe 


where E pe ,v pe and T are the Young’s modulus, Poisson ratio and thickness of the patches, 
respectively (it is assumed that the inner and outer patches have the same material properties). 
With <l;u , V\ and Vj denoting the piezoceramic strain constant and voltages into the two 
patches, the terms c pei = ty-Vj and e pe2 = ^ ly provide a relationship between the applied 
voltages and the resulting in-plane mechanical strains. 

From (4) and (5), it can be seen that the resultants due to the activation of the patches 
depend on the material and geometric properties of the patches, the radius of the shell, and 
the voltage being applied to the patches. Once determined, these resultants are substituted 
into (1) as the external loads on the shell. If multiple patch pairs are present, the resultants 
for each pair are incorporated in the shell equations in a similar manner. 


2.4. Normal loads and coupling conditions 


It can be seen in (1) that the shell equations also contain the term q n which represents the 
normal load on the shell. In writing the equations in this form, we have made the assumption 
that the only loads oil the shell are the contributions due to the patches, and the normal 
load which is due to the exterior noise field / as well as backpressure generated by the interior 
acoustic field. This assumption was made merely to simplify the discussion, and more complex 
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acoustic or mechanical loads on the shell can be incorporated as additional external forces and 
moments in the shell equations. Because q n , in this case, consists of contributions from the 
exterior and interior acoustic fields, the first coupling relation is 

q n {t,9,x) = f{t,6,x) - {t,w{t,6,x),v(t,0,x),u{t,6,x)) 

and is in general nonlinear since the backpressure onto the shell occurs at the shell s surface. 
The second coupling condition is the velocity constraint 

(t,w(t,6,x),v{t,Q,x),u(t,0,x)) = w t (t,6,x) 
or 

which simply states that the shell is impermeable to air. Note that this second constraint also 
provides a boundary condition for the acoustic response. 

2.5. Fully coupled system: strong form 

Under the assumption of small displacements which is inherent in the linear shell theories, 
the variable domain ft(t) is approximated by the fixed domain ft, and the general nonlinear 
coupling conditions are approximated by their linear approximations. Let the boundary Fo 
denote the shell in its unperturbed state. For s pairs of piezoceramic. patches, this then yields 
the linear model 

(j) tl = c 2 A(j> , (r, 9, x) G ft , t > 0 , 

V0-n = 0 , (r,0,x) € T ,< > 0 
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Recall that the internal shell moments and forces are summarized in (3) while the expressions 
for the external loads generated by the i th pair of piezoceramic patches are given in (4). The 
notation and [5i,2(0)]i denotes the indicator functions (see (2)) which are centered 

over the i tl1 pair of patches. Note that this representation admits different voltages and 
geometries for the ,s patches, thus increasing the flexibility of the model for control purposes. 

We also reiterate that the only damping in the system is the internal Kelvin-Voigt damping 
in the shell although medium damping inside the cavity can be added if cavity dimensions 
become significant. Without medium damping, however, the system is only weakly damped, 
and this must be considered when discussing theoretical convergence and well-posedness results 
for the problem. 


3 Weak Form of the System Equations 

As can be seen in the model (6), the use of the strong form of the system equations leads to 
first and second derivatives of both the internal and external moment and force resultants. 
This leads to difficulties both in approximating the behavior of the system and in solving the 
control problem. The first problem results from the presence and differing material properties 
of the piezoceramic. patches. As noted in [11] where piezoceramic patches are bonded to 
a beam, the material parameters of the combined structure must be modeled as piecewise 
constants in order to accurately match structural frequencies. Hence the parameters p, E , // 
and Cr> are expressed in terms of a Heaviside basis with the edges of the patches defining the 
support of the functions. This leads to difficulties in the strong form since it necessitates the 
differentiation of discontinuous material parameters. A similar problem arises when including 
the moments generated by the patches. As seen in (4), the support of the contributions is 
given in terms of Heaviside functions which implies that the use of the strong form yields an 
unbounded control operator (it involves the differentiation of the Heaviside functions as well 
as the Dirac delta). To avoid these difficulties, it is advantageous to formulate the problem in 
weak or variational form (the use of the variational form also admits the use of basis functions 
having less smoothness than those used when approximating the solution to the strong form 
of the equations). 

3.1. Weak formulation 

The second-order state for the problem is taken to be y = ( <j>,u,v,w ) in the Hilbert space 
H = Z 2 (H) x L' 2 (T 0 ) x L 2 {r 0 ) x L 2 (r 0 ). The choice of the space Z 2 ( Cl), defined as the quotient 
of L 2 (fi) over the constant functions, results from the fact that the potentials are determined 
only up to a constant. 

To provide a class of functions which are considered when defining a variational form of 
the problem, we also define the Hilbert space V = H'(Q.) x //^(To) x Hq(T 0 ) x Hq(T 0 ) where 
// 1 (0) is the quotient space of // 1 (Ti) over the constant functions. The subscript 0 in the 
remaining components of the product space denotes the subset of functions in the traditional 
Sobolev spaces which satisfy the essential boundary conditions u = v = w = |^ = 0 at 
x = 0, f. 
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A complete discussion concerning the derivation of the a variational formulation of the 
shell equations from energy principles can be found in [8]. For our purposes here, we simply 
note that integration in combination with the use of Green’s theorem yields the second-order 
variational system 
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for all ((,rji,T} 2 ,r} 3 ) € V (note that du = r dr d9 dx and d'y = RdOdx). The complex L 2 inner 
products are used in anticipation of the possibility of employing complex Fourier expansions 
in 9 when developing suitable approximation schemes for the problem. 

We note that, as is usual, in this variational form, the derivatives have been transferred 
from the plate and patch moments onto the test functions. This eliminates the problem of 
having to approximate the derivatives of the Heaviside function and the Dirac delta which 
arises in the case of the strong form of the equations. 

The system (7) can be formally approximated by replacing the state variables by their 
finite dimensional approximations and constructing the resulting matrix system. In order to 
discuss convergence results for the approximation, parameter estimation, and control prob- 
lems, however, it is advantageous to pose the problem in terms of sesquilinear forms and the 
bounded operators which they define. 


3.2. Abstract weak formulation 


Before proceeding with the abstract formulation, it is necessary to describe the inner 
products for the spaces H and V since these will be needed when determining the continuity 
of various operators. From energy considerations, it is appropriate to use the inner products 
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where $ = (<j>,u,v,w) and 'P = (£, rq , r? 2 , ?? 3 )- We remark that the inner product for the state 
space H contains terms which arise when considering the kinetic energy for the shell while 
the shell contributions in the V inner product are motivated by the form of the strain energy 
for a thin cylindrical shell (see [8]). Also, we note that with this choice of spaces and inner 
products, we can form the C.elfand triple V H ~ H m «-+ V* with pivot space H (see [27] 
for a complete discussion of these ideas). 

To define appropriate sesquilinear forms <r, : V x V — * (C , i = 1,2, we group the stiffness 
and wave contributions separately from damping and coupling terms, thus leading to the 
definitions 
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To account for the control contributions, we let U denote the Hilbert space containing the 
control inputs, and we define the control operator B <E C(U, V *) by 


<fl», ♦ .J t\ [(AT.U + J K*U £ - KW.U jg ~ K«.U W 


dm 
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for vp G V, where (-,-)v*,v is ^ ie usual duality pairing. Finally, by letting F = (°’°> 0 ’ i)’ 
we can write the system in the abstract variational form 

+ (T2(yt{t),ty) + = {Bu(t) + F t ^) v , y . ( 10 ) 

We reiterate that the state for this second-order system is y{t.) = (</>(/., •, •, • ),u(t , •, •), t »(/., •, •), 
«;(<,-,•)) in Vw ff. 

In order to write the system in terms of associated bounded operators, we first note 
that <T] and cr 2 are bounded (there exist c.\ and c, 2 such that |(Ti($,'P)| < c, |3>|v'|'P|v' and 
|cr 2 (<P, vp)| < c 2 1 ^ | v | 'P | v ) , and hence we can define operators Ai,A 2 G C{V,V*) by 

(A&,V) V .y = *,•($,#) 

for i = 1,2. This then yields the abstract second-order system 

yu(t) + A 2 y t (t) + A x y(t) = Bu(t) + F 


in V*. 

To apply infinite dimensional control results for periodic forcing functions to this problem, 
it is advantageous to write the system in first-order form. This is accomplished by defining 
the product spaces FL = V x H and V = V x V with the norms 

|(*, <1012, = |*|J, + 1*1?, 

and 

\(*,V)\v = Mv + I'Plv • 

We point out that V t — * Ft ~ FL* t -* V* again forms a Gelfand triple. 

The sesquilinear form a : V x V — ♦ (L' is then defined by 

a(Q, X ) = <r((T, A), (<J>, )) = - (A, <S>) V + oi(T, ) + «t 2 (A, t) 
where \ anc ^ © = (T, A). 

For the state 3^(0 = (y(0>2/*(0) we can subsequently write the system in the first- 

order variational form 


0 ? t(0>x>v*,v + ff 0W»x) = (Sii(f) + T(f),x)v,,v ( ll ) 

where jF(t) = (0, F(i)) and Bu(t) — (0, Bu(t)). As usual, the relation (11) must hold for all 
X G V- Finally, the weak form (11) is formally equivalent to the system 

W) = Ay(t) + Bu(t) + f(t) ( 12 ) 


in Ft where 


dom A = {0 = (T, A) € Fi : A € V, AiT + A 2 A G H} 


A = 


0 

—A\ 


l 

-a 2 


(id) 


ll 


Hence we see that the system corresponding to the original model can be written in various 
forms whose utility depend on the applications of interest. For approximating the dynamics of 
the system, the formulation (7), or equivalently (10) or (1 1), is useful and ultimately leads to a 
first-order matrix system when considering the resulting finite dimensional problem. Moreover, 
as will be shown in the next section, the weak form (10) provides a setting in which one can 
determine well-posedness results for the problem. Finally, the first-order infinite dimensional 
formulation (12) is an abstract Cauchy equation which facilitates the application of infinite 
dimensional LQR results to the problem of controlling interior acoustic pressure levels. 


4 Well-Posedness of the System Model 

In this section, the system model described in strong form in Section 2 and weak form in 
Section 3 is shown to be well-posed. The strategy used to do this can be summarized as follows. 
By employing the properties of the sesquilinear forms u\ and < 72 , it is demonstrated that the 
operator A of (13) generates a Co-semigroup S(t ) on 7i. Due to the fact that the control 
and forcing terms lie in V* rather than H, however, one is prevented from simply applying 
a variation of parameters approach to define mild solutions as is often done in problems of 
this type. Instead, a construction due to Haraux [17] is used to extend the semigroup 5(f) 
to the space W* = [dom .4*]* where an extension of the variation of parameters technique 
for describing a mild solution to the problem can be employed. By formulating the problem 
in this generalized manner, the well-posedness of the system solution can be established in a 
form that can be used in parameter estimation and control problems. 

4.1. Generation of the Co-semigroup 5(f) on H 

Before proving that the operator A generates a Co-semigroup on TC, we recall several 
definitions that will be needed (we note that these definitions vary slightly among authors, 
e.g. see [24, 27], and that the definitions given here will be used in the discussion which 
follows). 

We say that the sesquilinear form a : V x V —> C is V -elliptic if there exists a constant 
c > 0 such that 

Re a(4>, 4>) > c|$|y for all $ E V . 

The sesquilinear form a : V x V — * C is called H -semidliptic if there exists a constant b > 0 
such that 

Re <r($, 4>) > b\Q\ 2 „ for all $ e V . 

Finally, we say that the sesquilinear form a is symmetric if 

it( 4>, '!') = cr('P, 4>) for all $, 'I' £ V . 

The following theorem can be used to establish that the operator A given in (13) generates 
a Co-semigroup on 'H. The proof depends on the Lumer-Philips Theorem and arguments can 
be found iti [1] and pages 82-84 of [3|. 
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Theorem 1 . Let V t -+ H ~ H* <— » V* and suppose that <7] and <x 2 of (10) satisfy the 
properties: a\ is V-elliptic, continuous, symmetric, and a 2 is continuous, 77-semielliptic. 
Then A defined in (13) generates a Co-semigroup on H. 


To see that our system (10) satisfies the hypotheses of this theorem, we first consider the 
sesquilinear form <7](<I>, ty) = (<I>, ty) v given in (8). The boundedness of oq (which results from 
Schwarz’s inequality for inner products in conjunction with equivalence results for various 
Sobolev norms) was noted in the last section. The R-ellipticity and symmetry of <j\ follow 
directly from the definition of the sesquilinear form as the V norm. We note that the symmetry 
of u\ depends on the symmetry of the Donnell-Mushtari shell operator, and while some other 
shell theories such as that of Byrne, Fliigge and Lur’ye also provide symmetric operators 
and resulting sesquilinear forms, others such as the Love-Timoshenko theory do not yield 
symmetric operators (see [19]) and hence would not fall directly into our framework here. 
From a physical perspective, the symmetry of the shell operator and resulting sesquilinear 
form guarantees real vibration frequencies. This symmetry also guarantees that the shell 
model satisfies the Maxwell-Betti Reciprocity Theorem which essentially states that for a 
linearly elastic body, the general displacement at a point m resulting from a load at the point 
n is equal in magnitude to the displacement at n resulting from an equal magnitude load at 
m (see for example, [23]). 

We next turn to the sesquilinear form cr 2 of (9) and recall that the boundedness of a 2 was 
also noted in the last section. This follows from the definition of the o 2 and the equivalence 
of various norms which arise from the components of the V norm. The //-semiellipticity of 
c t 2 follows directly with the choice b — 0. 

The results from Theorem 1 guarantee that the operator A given in (13) generates a Co- 
semigroup S(t) on the state space 7i. Moreover, the semigroup satisfies the exponential bound 
|«S(f)| < e wt for t > 0 (where in fact, u> = 0 due to the fact that A is dissipative as shown in 

[3D- 

In the case of a bounded (in Ti) control input operator, the usual procedure is to use a 
variation of parameters approach to define a mild solution to the system. As noted in the 
last section, however, the control input B € £(£/, V*) defines the product space control term 
Bu(t) = (0,F?u(f)) 6 {0} x V* C V x V* = V*. Since Bu[i. ) lies in V* rather than in 7Y, the 
usual variation of parameters ideas are not feasible. We are therefore motivated to extend the 
semigroup S(t ) on H to a semigroup S(t) on a larger space W* D {0} x V* so that the mild 
solution 

w)=aw + /;#- S )( B u(s) ° +fw )* cm) 

is well-defined for Bu- |- F G L 2 (( 0, T ), V*) (in the following development, it will be shown that 
the space W* is given by W* = [dom^4*]*). This then provides a setting in which to guarantee 
the well-posedness (including continuous dependence on initial data and nonhomogeneous 
input terms) of the solution. 
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4.2. Description of A* and dom A* 


Before discussing the extension of S(t) from 7i to W*, it is useful to describe A * and 
dom A*. We first recall from (13) that A and dom A are given by 

dom A = {0 = (T, A) € H : A € V, A,T + A 2 A G H} 

• 40= ( -A t T-A 2 A ) ' 

From the usual definition of the adjoint A *, we now want to find dom A* and A* satisfying 
dom A* = {x = ($, ’f) € H|0 (A0)(x) is continuous on dom A } 

and 

M0,x)„ = (0,Xx)„ 

for all 0 = (T, A) € dom A and x = ($, € dom A*. It follows directly from the definition 

of the adjoint operator that x € dom A* if and only if there exists F = (71,72) € H = V x H 
such that 

Me,x)„ = (e,n„ 

for all 0 G dom A. The expansion of this relation yields the condition 

(A,<h) v + (Mi? - AjA), *)„ = (T, 7 i)v + (A, 72)// (15) 

for all 0 = (T, A) € dom A and x = € dom A*. By noting that (15) must hold for 

A = 0, it follows that 

(-A 1 T,'P) w = (T,7i)k = (AiT,7i)h 

for all T € dom A v The second inequality results from the observation that (T,7 i) k = 
o'! ( T , 7 ! ) = (v 4 iT, 7 ,) u . y = (A 1 T,7 1 ) h for 0 = (T,0) G dom A. It follows immediately that 
(A]T, — — 71 ) H = 0 for all T € dom Ai, which implies that 

7i = -V 

due to the V-ellipticity of o-j. Moreover, 71 € V implies that ^ £ V. 

When 7 j = is substituted into (15), the relationship 

{A] A, <f>) H + (—A^A, ty) H = (A,~f 2 ) H 

results for all A G V. From the last equality, it follows that the mapping A 1— > (A\A,$) H + 
( — A 2 A, ty) fj is continous in H. This in turn implies that — A^ G H and that 

(A, Aj4> — AW h — (A, j 2 ) h 

for all A G V. Hence we have 

72 = a:<d - a;v . 
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Thus A* and dom A* are given by 


dom A* = {X = ($, 9) € H\V € V, A** - e H } 


A*X = 




(16) 


A*4> - A* 2 ty ) 

Finally, we note that because £ V for any x € dom A *, it follows immediately that 

dom A* C V x V = V . 


4.3. Extension of <S(/.) to [dom A*]* 


With A* and dom A * described as above, we now consider the extension of the semigroup 
S{t ) generated by .4 from 77 to a larger space W* which is defined below. This will be 
accomplished by employing the extrapolation space techniques discussed by Haraux in [17]. 
Other, more detailed examples using these techniques can be found in [5, 6]. 

The space W = [dom A*] is taken to be dom A* with the inner product 


<*, *} w = <( A o - .4*)$, ( A 0 - 


for some arbitrary but fixed A 0 with A 0 > u> (recall that the original solution semigroup satisfies 
the bound |<S(f)| < e wt ). As proven in [6], the resulting W norm is equivalent to the graph 
norm corresponding to A*\ that is, there exist constants c i and c- 2 such that for any $ £ W, 
we have 


|(A 0 -^)<k| w < Cl (|4»| w + |X<&| w ) 


(17) 


|*|« + \A*<i>\ n < ci |(A 0 - A*)Q\ n . 


As discussed in [5, 6], the space W defined in this manner and with this norm is densely and 
continuously embedded in 77. Hence we can formulate the Gelfand triple W > 77 ~ 77* < — * VV* 
with pivot space 77. Moreover, the dual space W* is isomorphic to the completion of 77 with 
respect to the norm |ft|w* = |(Ao — A)~'h\n. 

To facilitate the arguments for extending the operator A from dom A C 77 to all of 77, we 
define the sesquilinear form a : 77 x W — * € by 


cr(0,x)s(eM*X> w 

for all 0 £ 77, x € W. Due to the equivalence of the W norm and the graph norm corre- 
sponding to A *, it follows that 

|<t( 0, x)l ^ l©|w|-4*x|w < c l©lw|\|vv • 


This then implies that for each 0 € 77, the mapping \ i-» <r(0, \) is in W*. Hence we can 
define AQ £ W* by (.40) (x) = <7(0, \) or equivalently, 

00) (x) = (0M-x>„ 
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for all 0 G W, X € W. From this definition, it follows that .4 G £(W, W*) with dom 4 = H. 
Moreover, if 0 G dom A, it follows from the Riesz representation theorem that 

(•40) (x) = (e,4'x)« = M0,x)k = (*40)(x) 

for all x € dom A *. Thus we see that >1 is an extension of the original operator A from 
dom A C 7X to all of 7i . 

Having extended the operator to the full state space, our final objective in this subsection 
is to observe that A is the infinitesimal generator of a Co-semigroup on W*. The proof of 
this result is given in [6], and we simply state the conclusion in the following theorem while 
directing the reader to that paper for details. 

Theorem 2. The operator A is the infinitesimal generator of a Co-semigroup S(t) on W* 
which is an extension of S(t) from TL to W*. 


Having described how the original semigroup S(t ) can be extended to a semigroup S(t) on 
W* = [dom*4*]*, we now want to show that the space {0} x V* which contains the control and 
force input terms is itself contained in W*. It is tempting at this point to recall that we have 
the Gelfand triples V ^ H ~ ^ V* and W ^ H ~ W W* as well as the embedding 

w (Z V, and from this assume that the desired result is obvious. We point out, however, that 
while W C V, the embedding is not in general continuous which implies that W c — > V can not 
be assumed to be true (if A* is V-elliptic for example, one can obtain the desired continuity 
of the embedding; however, this is not the case in this problem). Hence more care must be 
taken in order to argue that W* D {0} x V* 

To this end, let h = (0, A*) be an arbitrary functional in {0} x V*. In order to show that 
h E W*, we need to show that x = ($? > Mx) is bounded on W. Now, it follows that 

\h( X )\ = |A*(^)| < |A*|v* \ty\v < &]$|y . (18) 

Moreover, from the definition (16) of A* , we see that A*x — ( — — ^ 2 ^) which implies 
that 

mv < \A*x\ n < \A*x\ n + \x\n < c 2 \x\w • (19) 

(The latter inequality results from the equivalence of the W norm and the graph norm cor- 
responding to A * as described in (17),) From (18) and (19) it follows that there exists a 
constant c such that 

IMx)l < c \x\w 

therefore showing that h G W*. We have thus established the following result which is essential 
in developing expression (14). 

Lemma 1. Under the assumption of Theorem 1, we have {0} x V* C W* = [dom A*]*. 

In summary, we are now able to show the following. The system (11) is well-posed in the 
sense that (14) is a well-defined entity that can be taken as the mild solution to (11). More- 
over, because (14) is well-defined, the usual theorems for continuous dependence with respect 
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to initial data and nonhomogeneous terms (control and external forces) follow immediately. 
This facilitates the discussion of approximation and control ideas for these problems (see for 
example, [3, 4]). 

In the event that a 2 of (10) is also V-elliptic., one can argue (see [5] for details) existence of 
a weak or variational solution to the system (10) which agrees with the mild solution obtained 
from (14). In the example of Section 3, the form a 2 is not V-elliptic; nevertheless, one can 
turn to standard results from the theory of Co-semigroups, [22, 26], to obtain existence and 
uniqueness of solutions for the system under consideration in this paper. For example, we 
have the following result. 

Theorem 3. Consider the system represented by (10), (12) or (14) and suppose that the map- 
pings t. u{t) and t F(t) from [0, T] to 1R 1 and V*, respectively, are Lipschitz 
continuous. Then for each To £ BC = dom A., we have that (12) taken with T(0) = To 
has a unique strong solution given by (14). 

This theorem is readily established by appealing to Corollary 2.11, p. 109 of [22], in the 
context of our arguments above. Letting G(t) = (0, Bu[t) + F(t)) T and considering S{t) on 
the reflexive Banach space X = W*, we can argue that G is Lipschitz continuous under the 
hypotheses of the theorem. In particular, the inequalities (18) and (19) imply 

|G(<)|w* 5: c 2 \Bu{i.) T F(t.)\v , 

so that 

\G(t ) - G(s) |w* < c 2 {\B\ m v‘) |u(<) - u(. s)| + \F(t) - F(s)\v} • 

It follows that (14) provides the strong solution (i.e., differentiable a.e. in the W* sense - 
see [22]) to (12) interpreted in the W* sense. 


5 Conclusion 

In this paper, a modeling set of partial differential equations describing the dynamics of a 
structural acoustics problem was presented, and a mathematical framework amenable to the. 
development of approximation schemes for forward simulations, parameter estimation, and the 
application of PDE-based feedback control strategies was developed. The structure consists 
of a thin cylindrical shell with hard caps at the ends, and this component of the system was 
modeled by the Donnell-Mushtari shell equations. A constitutive law postulating that stress 
is proportional to a linear combination of strain and strain rate was assumed which yields a 
Kelvin- Voigt type of damping in the shell. The structural dynamics were then coupled to the 
interior acoustic field through pressure and momentum balance conditions. Finally, this model 
includes the contributions due to the activation of piezoceramic patches which are bonded to 
the shell and ultimately will be used to control the interior acoustic pressure levels. 

To accommodate the presence of the piezoceramic patches and differing material properties 
of these patches, the material parameters for the combined structure should be taken to be 
piecewise constants. This leads, however, to difficulties in the strong form of the system 
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equations since these parameters are contained in moment and force resultants which are 
differentiated when forming the equations of motion for the shell. Moreover, the discontinuities 
introduced by the patches lead to an unbounded control operator since it involves derivatives 
of the Heaviside function and Dirac delta. To avoid these problems, the weak form of the 
system equations was also presented. 

In the weak form, the derivatives appear on the test functions instead of on the moments, 
thus eliminating the difficulties associated with the discontinuous parameters and patch inputs. 
The weak form is also advantageous for many approximation schemes since it reduces the 
smoothness requirements for the basis elements. Finally, when the weak form is posed in the 
context of sesquilinear forms, convergence and well-posedness issues can be considered. 

The first step in proving the well-posedness of the system model was to verify that the 
first order system operator A : dom A C Tt — * 7i generated a Co-semigroup S(t) on Ti. This 
semigroup was then extended via the extrapolation techniques described by Haraux to a larger 
space W* = [dom A*]* so as to be compatible with the force and control inputs which lie in 
{0} x V *. Finally, the mild solution in terms of the extended semigroup was formulated and 
existence and uniqueness results were obtained. Hence the framework and model presented 
here is mathematically well-posed as well as amenable to the development of computational 
strategies. 
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